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The value of the dynamic critical exponent z is studied for two-dimensional superconducting, 
superfluid, and Josephson Junction array systems in zero magnetic field via the Fisher-Fisher-Huse 
dynamic scaling. We find z ~ 5.6 ± 0.3, a relatively large value indicative of non-diffusive dynamics. 
Universality of the scaling function is tested and confirmed for the thinnest samples. We discuss the 
validity of the dynamic scaling analysis as well as the previous studies of the Kosterlitz-Thouless- 
Berezinskii transition in these systems, the results of which seem to be consistent with simple 
diffusion [z = 2). Further studies are discussed and encouraged. 



I. INTRODUCTION 

The dynamics of two-dimensional (2D) zero-field sys- 
tems haxe-been studied continually over the last few 
decades JjlI usually Jel the context of the Kosterlitz- 
Thouless-BerezinskiiUna (KTB) transition. The dynamic 
critical exponent z characterizes the critical behavior of 
the dynamics of a transition. Despite the ongoing stud- 
ies of KTB dynamics, the value of z is usually not ques- 
tioned to be anything but the value that describes sim- 
ple diffusion: z = 2. This paper, through an analysis 
of various transport data sets from systems including 
superconductors (SC's), superfluids (SF's), and Joseph- 
son Junction Arrays ( J JAIs) using the dynamic scaling of 
Fisher, Fisher and Huse,tL9 presents ample evidence that 
the value of z in these systems is much higher: z ~ 5.6. 
The purpose of this paper is to convince the reader that 
despite the many previous reports consistent with z = 2, 
the question of the value of z is still an open one. Per- 
haps the single element that distinguishes this analysis 
from past analyses is that z is not taken a priori to have 
a value 2. In the dynamic scaling analysis, varying z to 
optimize the description of the data is standard. 

This paper is organized as follows. We begin in this 
section with a general description of the KTB transition, 
previous scaling attempts of I-V data, and previous find- 
ings of anomalous vortex diffusion. In Section ||, we will 
discuss the various length scales whose competition re- 
sults in the interesting critical behavior of this transition. 
In Section |nil we present scaling results on SC's, JJA's, 



The KTB transition is driven by the unbinding of vor- 
tex pairs. Below the transition temperature Tktb vor- 
tices are thermally induced and can only be excited in 
pairs that have a finite energy and not as "free" vortices, 
which have an infinite energy in an infinite system. As 
the temperature is increased, the number and size of the 
vortex pairs increase and these pairs screen one another's 
interactions. At the transition temperature, the vortex 
pairs start to unbind and free vortices are formed. It is 
the largest pairs that unbind first, leaving a finite density 
of smaller pairs above the transition temperature. The 
size of the largest pairs decrease as one goes further above 
Tktb- As we review below, for T > Tktb, the free vor- 
tices result in an ohmic I-V curve at low currents while 
the smaller pairs result in a non-linear I-V relationship 
at larger currents. 

Scaling techniques applied to zero-field I-V data from 
SC's or JJA's have been reported in the past but with- 
out finding the results reported here. In Ref. [ll], Wolf 
et al. scaled their I-V data from granular 30A thick 
films. There, however, the dynamic universality class 
was not explicitly studied— and no value of z was deter- 
mined. D. C. Harris et al.tn reported z — 2 in a dynamic 
scaling analysis of Josephson junction array I-V data, al- 
though it appears that they did no t allow z to vasy. We 
address that data in Section [II B. L. Miu et aZJl3 mcn- 



and SF's and also check for universality of the scaled 
data. We discuss the validity of the dynamic scaling ap- 
proach as well as that of the conventional approach in 
Section IV. We summarize the paper in Section 0. 



tion peripherally a dynamic scaling analysis on BSSCO 
(2223) in which it is found that z ~ 4.5 (assuming a 2D 
system). However, those results were not pursued fur- 
ther those authors. In the case of superfluids, Brada 
et al.[3 have performed a finite size scaling analysis (in 
contrast to a dynamic scaling analysis) of their frequency 
dependent superfluid density and dissipation data from 
helium films. They however did not include a study of 
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the value of z in their work. 

There has been experimental work that popts to- 
ward anomalous vortex diffusion. Theron et alS3 used 
impedance measurements on weakly frustrated JJA's to 
provide evidence of non-conventional vortex dynamics. 
Unfortunately, this work did noL-probe the value of..*- 
There has also been theoretical^ and simulationaBEZl 
work that indicates anomalous vortex diffusion have been 
made, but none conclude that z > 2 at T — Tktb- In 
this work, in order to limit the scope of this already long 
paper, we will not discuss the theoretical and simula- 
tional work, but concentrate only on experimental data. 

In comparing the "conventional" approach with the 
dynamic scaling approach, we will make use of many of 
the formulas derived from the former. However, because 
the dynamic scaling approach indicates that the dynam- 
ics are non-diffusive, the assumptions used to derive the 
formulas for the dynamics in the conventional approach 
may be incorrect. For that reason, we will use those for- 
mulas only to address the validity of the conventional 
approach and not the validity of the dynamic scaling. 

It is useful for us to clarify our use of the phrase "con- 
ventional approach." Because the dynamic scaling has 
only been used sparingly in analyses of the type of data 
that we look at here, we view this approach as non- 
conventional. We will therefore refer to any other ap- 
proach to analyzing this data as "conventional" or tradi- 
tional. 

The value of z that is found here runs contrary to some 
of the conclusions of previous studies. However, we stress 
that our results are not inconsistent with all of the results 
in the literature. In particular, the dynamic scaling re- 
sults presented here do not contradict (or even pertain 
to) the conventional findings for static behavior of the 
KTB transition. Further, our results differ from only 
some of the experimental st udies of the dynamics. This 
will be discussed in Section IV B but we mention some 
of those now. If the measurement does not involve z 
or if the measurement does involve z but can measure 
only the product bz (where b is a material-dependent 
constant that enters through the correlation length [see 
Eq. (Q)] ) , then our findings do not contradict those mea- 
surements. An example of the first type of measurement 
are "static" kinetic inductance measurements. Ea Exam- 
ples of the second class of measurement are the resis- 
tance, noise spectrum-4peasurements,E3 and helium tor- 
sion measurements. E2IEII 



the KTB transition and critical behavior as well as the 
approaches used to study them is needed. In this section, 
we will review this background, directing our discussion 
primarily at superconducting systems. (For discussion on 
JJA's or superfiuids, see the reviews listed below.) We 
will discuss the various relevant length scales, the crite- 
ria for a phase transition, and the approaches one can 
take to study the KTB dynamic critical behavior. We 
don't intend our review to be comprehensive and j^fer 
the reader to any of the many excellent reviews.LJO 



A. Length Scales 

The competition of the length scalesH in the system 
determines the critical behavior of the KTB transition. 
For this reason it is important to review each of them. 
One can subdivide the length scales of the system into 
two categories: intrinsic and extrinsic. By extrinsic, we 
mean those length scales that are determined by an ap- 
plied current or magnetic field. 

Intrinsic 

The intrinsic length scales include the vortex correla- 
tion length £(T), the 2D penetration depth X2D — 2 A 2 /d 
(where A is the London penetration depth and d is the 
sample thickness), and sample size. 

There are three important aspects of the correlation 
length that we discuss here. The first is its distinc- 
tive temperature dependence for temperatures above the 
transition temperature Tktb'm 



£+(T) <x e X p[Vb/{T/T KTB - 1)] 



(1) 



where b is a non-universal constant. This unique temper- 
ature dependence is in contrast to the common power-law 
dependence one finds, for example, in Ginzburg-Landau 
theory. 

The second aspect is the behavior of £_(T) below 
Tktb- Because the susceptibility below the transi- 
tion temperature is infinite, Kosterlitz originally defined 
£-(T) to be infinite. Based on the critical behavior of 
the dielectric constant, Ambegaokar et ril. (AHNS)EJ de- 
fined a finite diverging correlation length for T < Tktb- 
The two results do not contradict one another since they 
have different meaning. The AHNS correlation length for 
T < T^TB-pan be thought of as_the size of the largest vor- 
tex pairs.ES Ambegaokar et alu estimate that the £-(T) 
has a smaller magnitude than £ + (T): 



II. THE KTB TRANSITION IN 2D 
SUPERCONDUCTING SYSTEMS: 
BACKGROUND 

In this paper, we propose an interpretation of trans- 
port data on superconductors, superfiuids and Josephson 
Junction arrays that is very different than that which 
has been accepted for the last twenty years. In order to 
judge the two approaches, a thorough understanding of 



£_ (T) (x exp [VVMl - T/Tktb)] 



(2) 



In this paper, we will take £_ (T) to represent the size of 
the largest pairs. 

The third aspect of the correlation length that is im- 
portant in this paper is its behavior in an applied current 
/. The effect of an applied current is to unbind vortex 
pairs down to zero temperature and therefore to destroy 
the phase transition. As a result, the correlation length 
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no longer diverges for finite I at Tktb and has the fol- 
lowing behavior: 



£±(T,/)ocy/[J£±(T,/ = 0)/T], 



(3) 



where / is a non-singular function. 

The next two intrinsic length scales are associated with 
finite size effects and give a single vortex a finite "bare" 
energy. The first, which does not apply to superfluid 
Helium, is the 2D penetration depth X2D — 2\ 2 /d. At 
distances less than this length from the vortex core, the 
superfluid velocity goes as 1/r. Beyond, this length, the 
superfluid velocity decreases as f/r 2 .E3 In a "perfect" 
(d = 0) superconductor, A2D = 00 and the superconduc- 
tor behaves as a superfluid would. The second finite size 
length is sample size, which is the smaller of the sample 
width W and length L. (Typically W <C L.) The energy 
of a free vortex is, 



E FV = [q 2 /2] ln(L fs /£ ) + E c 



(4) 



where £0 is the size of the vortex core and the finite size 



Lf a = min[A 2 D, W]. 



(5) 



The vortex interaction strength q 2 and the core energy 
E c depend on the system (SC, SF, of JJA). A finite en- 
ergy for a single vortex means that there will be free 
vortices below the transition temperature, which in turn 
precludes a true phase transition as we will discuss in the 
next section. 

Extrinsic 

Three extrinsic length scales characterize the applica- 
tion of dc and ac electric and magnetic fields. The first is 
the length scale r c ,HlH3 which is the length scale probed 
by a dc applied current. For a superconductiattJilm, the 
energy of a vortex pair with separation R iscirESI 

E(R) = [Trn 2 s D h 2 /2m] ln(R/£ ) - irMRd/eA + 2E C (6) 

where the vortex interaction strength has been ex- 
pressed in terms of the superconducting parameters q 2 = 
Trn 2D h 2 /2m (n 2D — n s d is the areal superfluid density, 
n s is the superfluid density, and m is the mass of a 
free electron,) and A{= Wd) is the cross-sectional area 
through which the current I flows. For small separations 
(R < r c ), the 2D, logarithmic term dominates and the 
interaction is attractive. For R > r c , the linear, current- 
induced term dominates and so the interaction is repul- 
sive. As a result, the interaction energy peaks at 

r c = 4k B T KTB We/irhI (7) 

where we have written the interaction strength in 

^ 2D h 2 /2m = 



terms of the transition temperature,!! 
AksTxTB- The non-linear I-V relationship originates 
in thermally-activated hopping over this barrier at a rate 
r which depends on the value of the vortex pair energy 



E(R) at this separation: T oc exp[E(r c )/k B T]. Thcm- 
fore, the dc I-V curves probe length scales of 0(r c )Bt3 
In ac measurements (e.g. kinetic inductance) with cir- 
cular frequency u>, the probing length is the diffusion 



length, r u = 
sion constant 



/(x)) 1 / 2 , where D is the vortex diffu- 
This result is derived by analyzing 
the linear response of the dielectric constant. It should 
be pointed out that the two quantities, r c and r w , do not 
compete per se with the other length scale but rather 
indicate the length scale being probed. 

The final extrinsic length scale that we mention is 
due to an applied magnetic field and characterized by 
the average distance between field-induced vortices l B — 
(&0/B) 1 / 2 where <E>o is the superconducting flux quan- 
tum and B is the magnetic induction. The field-induced 
vortices are "free" and present at all temperatures, which 
precludes a true phase transition. 



B. Existence of the phase transition? 

Because of X2D, Kosterlitz and Thouless originally 
wrote that this [-critical behavior would— pot apply to 
superconductors.Q It was later realizedcSEII that in prac- 
tice, A2D can be larger than the system size and so su- 
perconducting films should not behave much differently 
than superfluid films. Whether A2D is larger than the 
system size or vice versa, there will be a finite density 
of free vortices below the transition temperature. The 
density of free vortices will beau 



T q 2 /2k B T 

n F cx LfJ 



(8) 



where q 2 here is the renormalized vortex interaction 
strength and, again, Lf s = min[A2£> , W] . Note that sys- 
tem size and X2D enter Eq. (||) in the same manner, im- 
plying that one cannot tune the ratio of X2D to W to 
observe the transition. This is contrary to a perception 
in the literatureEHTEj that the transition can be observed 
if A2D > W. This begs the question of why one should see 
critical behavior at all if there is no true phase transition 
in finite size systems. The answer of course is that critical 
behavior can be seen if there is a diverging length. This 
occurs in the KTB system, provided that £±(T) < Lf s . 



C. Approaches to studying KTB dynamics in 
superconductors 

A variety of approaches can be used to study KTB dy- 
namic behavior. Here we review the two main approaches 
used to study the dynamics of superconducting and JJA 
films and that therefore determine a value for the dy- 
namic critical exponent z. Brief mention of the methods 
used for su perfl uid hel ium systems will be made in Sec- 
There are other approaches for 



IIIC and IV B 2 



tions 

investigating KTB behavior that don't uniquely deter- 
mine a value for z: we will not review here but we will 
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discuss them in Section IV B in the context of previous 
evidence that z — 2. 

In this section, we will first review the conventional 
results and its derivation. We will then generalize these 
formula to a general value of z. The dynamic scaling 
analysis will then be introduced and finally, the connec- 
tions between the two approaches will be discussed. 



R oc n 



1. Conventional approach 



In the first, more "conventional" 
current-voltage (I-V) isotherms are-jp 
alyzed in terms of their 7 — > limit 

V/Iozl^- 1 . 



approach, the 
iasured and an- 



(9) 



The signature of a KTB transition is a jump from 
non-linear behavior below the transition temperature to 
ohmic above the Txtb'- 



V/I = R(T) oc exp[-2y/b/(T/T K TB - 1)] (10) 

where b is a non- universal constant. [Minnhage n0 has 
generalized this to take into account the underlying su- 
perfluid. See Eq. (p2|).] In particular, the exponent a(T) 
will decrease linearly with increasing temperature until 
the transition temperature is reached at which point a(T) 
will, in the 7^0 limit, jump from 3 to lEj because of 
the "universal jump" in the superfluid density.E£l Indeed, 
the condition a(T — Tktb) — 3 is commonly used to 
determine the transition temperature. 

Eq. (Q) is derived by determining i?(T, 7) above and 
below T KT B and usiagJZ = IR(T,I). This derivation 
is well documentedcjQ'Ej but we will highlight the key 
points. To find R(T,I), the density of free vortices n/ 
must be determined. 

Below the transition temperature, free vortices in the 
limit of a weak current are produced by thermal activa- 
tion over the barrier in E(r c ) as mentioned above [below 
Eq. (H)] . This is done using the kinetic equation for the 
rate of change in the number of free vortices: 



dnf/dt = T(T, 7) - 



(11) 



As mentioned above, T(T, I) is the rate at which vortex 
pairs are unbound [TjT, 7) oc exp{-7;(r c )/£; B T}.] The 
second term in Eq. (|ll| ) takes into account free vortices 
combining to form pairs. In steady state, 



n F =r(T,7) 1 / 2 cx7« 2 / 2 ^ T , 
Above the transition temperature, 



(12) 



(13) 



since £ + is the average distance between free vortices. 

The final step in datejfHjining Eqs. (p|) and (|l0|) is to 
relate R(T,I) to n/l 8 » 



(14) 



Substituting Eq. (0) and Eq. @ into Eq. (P) and 
using V = IR(T, I) onej-arrives at Eqs. (||) and (|10|) re- 
spectively. Kadin et al.c3 have made extensions of this 
work to finite current. Many workers in the field however 



take Eq^ 

Eq. 
formula] 



to be valid over wide ranges of 7. 
) is based on the Bardeen-Stephen. flux-flow 
(as stated by subsequent authors.E3o) For 



subsequent discussion, it will be useful to outline the 
derivation of this equation,E3 whose starting point is the 
electrodynamic Josephson relation, 



V = [h/2e]dA(f)/dt, 



(15) 



where A0 is the change in the phase of the supercon- 
ducting order parameter across the width of the sample. 
dAcj)/dt is proportional to the number of vortices that 
cross the width of the sample per unit time, 



\dA(j)/dt\ = 2nLn F \v D \ 



(16) 



where vd is the vortex drift velocity. Finally, one assumes 



I'd = [xirhl/eA, 



(17) 



where fi is the vortex mobility and e is the electron 
charge. Because the vortex mobility is taken to be local 
and therefore independent of tip, Eqs. (|l5|)-([l7|) result in 
Eq. @: R=V/I(xn F . 

The linear relationship between R(T) and nj in 
Eq. (|lj) presumes single vortex diffusion. If one were 
to allow more complicated dynamic or critical behavior, 
then Eq. ( |l4| ) must be modified. The general expression 
for R(T) in the critical region depends upon the dynamic 
exponent z: 



RocC 



(18) 



Of course the two expressions are identical for z = 2. 
Using Eq. r^F pl ace of R oc rif in the conventional 
derivationJiSEj'Ej one would find that a(T) will jump 
from z + 1 to 1 in the 7^0 limit and that 



R(T) oc expi-zy/b/iT/T, 



KTB 



!)]■ 



(19) 



We will discuss t he ram ifications of using z / 2 in 
Eq. (pi) in Section [TIC 3. 



2. Dynamic scaling approach 

Dynamical scaling is motivated by the observation of 
critical slowing down near a continuous transition. This 
phenomenon is well established for the KTB transition, 
and is marked by the divergence of the relaxation time 
scale, r. The dynamic scaling hypothesis asserts that 
critical slowing down is related crucially to the divergence 
of the static correlation length£3 t«( z . 
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In general, many types of dynamics may be associ- 
ated with a particular static universality class, and thesa 
should fall into distinct dynamical universality classes.E2l 
For two-dimensional systems, including SC's, SK's, and 
JJA's, the conventional KTB dynamical theoryoH is con- 
sistent with model A universal dynamics ,Ej and z = 2. 
For bulk superconductors however, the correct dynami- 
cal universality class is_presently unclear, and seems not 
to belong to model A.eII 

In light of this situation, it is important to test the 
conventional theory by performing a scaling analysis, in 
which the dynamical exponent, z, and the scaling func- 
tions are a priori unspecified, z is then determined by 
collapsing different data sets onto a single curve. Al- 
though conventional KTB dynamical theories do not 
support such a general approach, a scaling ansatz has 
recently been proposed by Fisher, Fisher, and HusetHl 
(FFH) for superconducting critical phenomena. This 
successful ansatz has been applied to a wide variety of 
systems and transitions, including simulations of KTB 
dynamics .E3 We expect that dynamical scaling will be 
appropriate in the same regimes where the static KTB 
theory is applicable. 

In the FFH theory, for 2D superconductors, the I-V 
curves should scale as 



v = ic z x±(WT), 



(20) 



where x+(-){ x ) is ^ ne scaling function for temperatures 
above (below) Tktb- The two important asymptotic be- 
haviors of x( x ) are um 2;^o X+i x ) — const, (ohmic limit), 
and linxr^oo x±i x ) ^ x * (critical isotherm). The uni- 
versal jump appears as the difference in (log-log) slopes 
between the two asymptotic limits of x+{ x ), as we will 
discuss below. 

It is convenient to rewrite the Eq. ( |20| ) as 




(21) 



where e±(x) = x/x± ( x )- The advantage of Eq. ( Bj ) 
over Eq. (20) is that one can better judge the scaling, 
because only the x-scale is stretched in Eq. (^l]). In 
Eq. (^0|) , both the x-scale and y-scale are stretched mak- 
ing it harder to judge a collapse of the scaled data. (Com- 
pare Fig. 4 of Ref. |l| with Figs. [fj-| here.) 



3. Connections between the two approaches 

While not apparent at first glance, the connections be- 
tween the dynamic scaling approach and the conventional 
approach become clear when one considers the following. 
First, one must keep in mind that t he conventional ap- 
proach, as described in Section II C 1 , is valid only in the 
limit / — > while the dynamic scaling approach is valid 
for finite currents. Secondly, the relationship of one to 



the other should be considered for an arbitrary value of z. 
[For the generalization of the conventional theory to an 
arbitrary value of z, see the discussion around Eqs. ( |lg| ) 
and ©]. 

Taking these considerations into account and looking 
at the asymptotic limit / — > of the dynamic scaling 
function, one finds that the two approaches are indeed 
compatible: (i) both theories predict that the critical 
isotherm (T = Tktb) should be a power-law V cx I z+1 
[i.e., a{TxTB) = z + 1], (ii) for T < Tktb, both theories 
agree that the voltage remains a power-law of the current 
(iii) for T > Tktb, both theories give a(T) — 1, with 
R{T) defined as in Eq. @. 



III. SCALING RESULTS 

In this section we will apply the scaling theory 
[Eq. (|2l]) l of Fisher, Fisher, and Huset3 to transport data 
from superconductors, Josephson Junction arrays and su- 
pcrfluids. The universality of the scaled data will then be 
checked. Preliminary results on both high-temperature 
superconductors (HTSC's) and low-temperature (con- 
ventional) superconductors (LTSC's) were presented in 
Ref. |^| where it was concluded that z ~ 5.7 ± 0.3. 

A preliminary estimate for the value of z indicates 
that z> 2 for these systems. The critical I-V isotherm 
(T = Tktb) is easily identified on a log- log plot by the 
fact that it is straight (i.e., ol{Tktb) is independent of 
current). The slope of this isotherm [ol{Tktb) = z + 1] 
gives an estimate for z. A visuaLcheck for this condition 
on the TV data-of Repaci et al.pj Vadlamannati et al.p3 
Matsuda et al.jtB, or any of the others clearly shows that, 
indeed, z ^> 2. With this initial evidence, we move on to 
a more rigorous scaling analysis of the data. 

To perform the scaling analysis, £± (T) in Eq. ([H]) must 
be specified. This can be done for superconductors and 
Josephson junction array by exploiting the ohmic limit 
of Eq. (H): R(T) cx £+ z - (For superfluids, the thermal 
conductance K is used in place of the resistance: K cx 
For we will assume that the vortex correlation 
length is symmetric about the transition (modulo some 
pre-factor) in this section. We will explore the validity 



of this assumption in Sec. [V A 1 



Note that in the dynamic scaling theory, there are no 
requirements for the temperature dependence of £±(T). 
In this work, we will assume that the KTB form, £(T) cx 

exp[y/b/ (T/Tktb — 1)], provides the most efficient pa- 
rameterization of the correlation length. It is through 
this assumption that the explicit connection with KTB 
theory is made. Any other temperature dependence for 
R{T) could be used to check the scaling collapse, which 
would leave the type of transition more ambiguous. 

We determine R(T) in two ways where possible. The 
first method is to extract it from the ohmic part of the 
I-V curves and the second method is to digitize the R(T) 
data. The two results are then compared to one another. 
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In the case of discrepancies, the R(T) determined from 
the I-V is used since thermal equilibrium is more likely in 
that case. Another advantage of using R(T) determined 
from the I-V is that one is assured that the R(T) is 
ohmic. The disadvantage of course is that fewer data 
points are available for this R(T). It should also be noted 
that we fit R(T) only over the temperature range over 
which we have I-V isotherms. 

The fitting parameters for Eq. ( pl| ) are z (universal), 
Tktb and b (non-universal). Three requirements must 
be fulfilled self-consistently in our scaling procedure: (i) 
V oc I z+1 , along the critical isotherm T — Tktb', 

(ii) R(T) oc £ -z , in the high temperature range; and 

(iii) scaling collapse of the I-V isotherms, according to 
Eq. (12]]) ■ Condition (i), which says that the I-V curves 
are straight on a log- log scale at Tktb, is used first to 
estimate a value of Tktb and z. That value of Tktb 
is then used in (ii) to fit the ohmic resistance data to 
obtain an expression for £± (T) . Finally, condition (iii) is 
checked. Because there may be a couple isotherms that 
appear to be straight on the log-log scale, this process is 
repeated, in the manner of Shaw et alrB, for the accept- 
able range of Tktb's to satisfy all three conditions. 

In some of the data sets that we examined, the I-V 
curves crossed over to an ohmic region at large I.E3 This 
behavior is not due to the critical behavior of the vortices 
but rather to a breakdown of the underlying superfhiid. 
For that reason, we have omitted such data from our 
analysis. 

Note that in the following, we have displayed all the 
scaling results (except for a few noted exceptions) with- 
out first making a judgment of the quality of the data. As 
a result, the quality of the scaling also varies. Nonethe- 
less, we stress that each scaling result displayed here has 
been optimized for the best collapse and not for a value 
of z in agreement with the other samples. This makes the 
result that all of the collapses occur for roughly the same 
value of z all the more striking. Furthermore, we believe 
that there is a strong correlation between the quality of 
the raw I-V data and the scaled data. But we leave this 
for the reader to judge. 



A. Superconducting films 

Fig. pi s hows the scaling of three separate data sets 
to Eq. (pl|), (previously reported in Ref. [l3] but plotted 
on a different scale here.) The first data seta (marked 
a in Fig. [I]) covers a temperature range [30.06 K:46.09 
K] and is from a YBCO/PrBa2Cu307_<5 multi-layer in 
which the YBCO layers have a thickness of 24 A and 
are reported to be nearly electrically isolated from one 
another by PrBa2Cu30y_5 (PBCO) barrier layers. The 
scaling procedure leads to the results Tktb — 32.0 K, 
b = 14.0, and z = 5.6 ± 0.3 where the resistance was fit 
over the range [43.5 K:47.0 K], (See Figure Jh Curve b 
in Fig. [l] is the scaled I-V data from Ref. 34 taken on 



a 12 A thick YBCO mono-layer and includes isotherms 
ranging in temperatures from 10K to 40K. The resistance 
was fit over the range 25K to 34. 5K yielding, along with 
the other two criteria, the following parameter values: 
Tktb = 17.6 K, b = 7.79, and z = 5.9 ± 0.3. The 
scaled data set denoted by c in Fig. |l| and covering the 
temperature range [2.6 K:3.4 K] corresponds to_a con- 
ventional, 100 A thick, superconducting sampleEEl (Hg- 
Xe alloy). The parameters which led to this collapse are 
Tktb = 3.04 K, b = 3.44, and z = 5.6 ± 0.3. The resis- 
tance was fit over the temperature range [3.3 K:3.4 K]. 
While this collapse is not as complete as those of the 
others in this figure, we emphasize that the best collapse 
was obtained for the reported value of z. 
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FIG. 1. The I-V curves from thin superconducting films 
scaled with Eq. @: (a) 24 A thick YBCO layars in a 
multi-layer structure from S. Vadlamannati ret a/.;ufl (b) 24 
A thick YBCO mono-layer from Renap et al.jtH and (c) 100 A 
thick Hg/Xe film from Kadin et alxS (Data sets (b) and (c) 
have been shifted arbitrarily.) The lower branch of these plots 
correspond to T > Tkt b ■ The limits of this branch are ohmic 
in the weak current limit to V oc I z+1 in the high current 
limit. It is these limits that represent the jump in the expo- 
nent a. 

A few features of the scaled data in Fig. [j] should be 
pointed out. First, the upper branch corresponds to tem- 
peratures below the transition temperature and the lower 
branch to temperatures above the transition tempera- 
ture. Secondly, an ohmic TV relationship here is repre- 
sented by a slope 1 on the log-log scale. One can see in 
each of the three curves in this figure (and the scaled data 
in the following figures,) the lower (T > Tktb) branch 
is ohmic at low values of the scaling variable x (typically 
low currents) and curves over and approaches a horizon- 
tal line as x is increased. The horizontal line corresponds 
to the I-V relation, V oc I z+1 . 
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FIG. 2. The scaled I-V curves of (a) Ref. |§ on a YBCO 
mono-layer; (b) R ef. k6 on a In/InO 150 A thick composite 
film; and (c) Ref. f|jTrom a 1000 A thick BSCCO crystal. 
(Data set (b) and (c) has been shifted arbitrarily.) The col- 
lapse of curves (a) and (b) in this figure are not as complete 
as that of curve (c) or curves (a) and (b) of Figure [j]. Yet, 
the collapse could not be improved using other values of z. 



scale in the same way as the thinner samples in Section 
IV A3 . 

of the Hebard and Fiory I-V data 
was not possible since the temper- 



A scaling anahisi; 
on a In/InO filmliac 
atures of their isotherms were not published. 



B. Josephson junction arrays 

I-V characteristics of Josephson junction arrays are ept. 
pected to be similar to that of superconducting films.Ej 
In this section, we apply Eq. (|T]) to data from two JJA 
systems. Etlj A primary difference for these systems is 
that their resistance is not described by Eq. ( [l9| ) be- 
cause the temperature is renormalized and depends upon 
the temperature-dependent critical current i c (T). In the 
data that we analyze below, we did not have access to 
i c {T) and so we could not determine if R(T) followed the 
KTB behavior for JJA's. Nonetheless, following the dis- 
cussion at the beginning of this Section (III), we will use 
the Eq. fll9[ ) to check for a collapse of the data according 
to Eq. (|21|). As a result, the value of b will not have the 



significance it had in Section III A 



For this paper, we have extended this analysis to two 
more 2D superconducting films. The collapse is shown in 
Fig. H The curve marked a in that figure is from a mono- 
layer of YBCO sancbajiched between two PBCO layers of 
different thicknessesc3 and covers isotherms ranging from 
16.32K to 41.31K. For the parameters, we find Tktb — 
18.3 K, 6 = 31.04, and z = 5.3 ± 0.5. The value of T KTB 
is similar to that of the YBCO mono-layer of Ref. |34| but 
the value of b is nearly 4 times larger. Curve b in Fig. || is 
from I-V data on 150A thick In/InO composite film from 
Ref. ^6|. The I-V isotherms cover a temperature range 
[3.010 K:3.182 K]. Here, the resistance data determined 
from the I-V curves covered only a limited range and did 
not match well the R(T) data from Fig. 3 of that paper. 
For these reasons, the scaling criteria (ii) was not fully 
met. Nevertheless, the scaling collapse occurred for a 
value of z near that of the other samples: Tktb = 2.97 K, 
b = 10.21, and z = 5.2 ± 0.5. _ 

We also applied this scaling analysiscS to a thicker 
(1500 ± 500 A) Bi 2 Sr 2 CaCu 2 08 + 5 crystal. For such a 
thick crystal of a layered material, one would expect a 
crossover to 3D behavior near the critical temperature 
and a failure of the 2D scaling and perhaps a breakdown 
in the scaling. We however did not see any breakdown in 
the scaling for this sample. Indeed, as shown in curve c 
of Fig. ^, we found a good collapse of the data with the 
2D scaling form with z = 5.6±0.3, T KTB = 78.87 K, and 
b = 0.20. (See also Fig. 2 of Ref. ||) A difference be- 
tween this scaled data and the others was found when we 
looke d at universality. We examine this issue in Section 
HID and discuss the reason why a thicker crystal may 



10"' 



3 



3. 10 




10"'" 

10"" 10"" 10" 10" 

l£/T (arb. units) 

FIG. 3. The scaled I-V curves from Josephson junction ar- 
rays: (a) YBCO/Ag weak links from Ref. g (b) Nb-Ag-Nb 
proximity-coupled junctions from Ref. |l§ and (c) Nb/Nb ar- 
rays of Ref. [48l (Data sets (b) and (c) has been shifted arbi- 
trarily.) 

Curve a in Fig. ||| is scaled I-V data from an high- 
temperature_superconducting Josephson junction array 
(YBCO/Ag)c3 where the parameters used were Tktb — 
74.3 K, z = 5.8 ± 0.4 and b = 0.72. The resistance was 
fit over the temperature range: [78K:90 K], One can see 
that the collapse is very good except for the isotherms 
furthest below Tktb- This breakdown could be because 
those isotherms are out of the critical region. An attempt 
to optimize the T -C Tktb collapse by letting the corre- 
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lation length be asymmetric in accordance with Eq 
was unsuccessful. 




The scaled data denoted by curve b in Fig. 
a Nb-Ag-Nb proximity-coupled junction array, 
used the "100%" data from Ref. ^0[) The parameters 
which produced the best fit are Tktb = 6.84 K, z = 
5.8 ± 0.3 and b = 0.32 and the resistance fit was over 
the temperature range: [7.3K:7.8 K]. The scaled I-V's 
covered the range: [6.9K:7.8 K] . One can see that he 6.9K 
data set which does not scale well for low currents. We 
believe that this is not a real effect since all of that non- 
collapsing data has a voltage V < 10 _9 V^ and does not 
parallel the behavior of the other data in that rang©. The 
data do scale well for V > 10~ 9 V. Harris et a/.l!3 have 
previously checked the dynamic scaling of this data and 
concluded that z = 2. We believe that they would have 
reached the same result as found here had they allowed z 
to vary to optimize the scaling and had they used Eq. ( |2l| ) 
instead of Eq. (£0 ) . 

Curve c in Fig. | is the-scaled I-V data from a Nb/Nb 
Josephson junction array.Ej The parameters used to opti- 
mize this were tktb = 0.51 K, z = 5.5 ±0.5 and b = 5.7, 
where r = ksT /[(h/2e)i c (T)]. A collapse of the scaled 
data could be obtained over a relatively large range of z 
and Tktb- We attribute this to the fact that the J-V's 
isotherm only covered 2-3 decades of voltage and that 
there were no isotherms for temperatures below Tktb- 
The resistance was fit over the range r = [fLS:1.3]. 

In a conventional analysis, Abraham et al.c3 reported a 
jump in the TV exponent a of 3 to 1 for PbBi/Cu arrays, 
a result which implies z = 2. As far as we could tell, that 
data was not published and so was not available for our 
analysis. We tried to apply Eq. (|2l]) to sample 6-18-81 of 
Ref. ^2] but no definitive conclusions were reached due to 
the fact that the I-V isotherms only covered 2-3 decades 
of voltage and a limited temperature range. 



C. Superfluid He films 

For superfluid 4 He films, the analog of electrical con- 
ductance and I-V curves are thermal conductance and 
Q-AT curves, where Q is power through the film and 
AT is the temperature gradient across the film. These 
measurements are done by injecting heat at one end of 
a thin superfluid film adsorbed on a surface (e.g. my- 
lar). Because of dissipation from vortex pairs, there is 
a temperature gradient across the film that is measured 
and is the analog of the voltage in the superconducting 
measurements. In reality, Q is not the heat through the 
film but rather represents 4 He mass flow from the cool 
end of the film to the warm end, which occurs to replen- 
ish the 4 He which evaporates from the warmer end at a 
quicker rate than from the other end. This is thoroughly 
discussed in Ref. |||. 

As we shall shortly see, these results do show that z ~ 
5.6. However, we are careful to point out that the results 



we are about to present can only be said to be consistent 
with such a value of z but cannot be taken to be evidence 
that z ~ 5.6 for two reasons. First, there is no reliable 
thermal conductance K data in the "ohmic" limit (i.e., 
T > Tktb zero power Q limit.) Furthermore, one has to 
account for the conductance of the gas K g surrounding 
the film which can only be approximated to within a 
factor of two. These two points allow us more freedom to 
obtain the optimal scaling. Nonetheless, the best scaling 
does yield a z which is similar to that of superconductors 
and Josephson junction arrays. 

The collapsed data in Fig. [| marked curve a is from 
Ref. [33]. Instead of varying the temperature, those au- 
thors varied film thickness d which in turn varies the 
transition temperature. The independent variable then 
is d and not T and so the correlation length depends 
upon \d — d c \ instead of |T — Tktb\ where d c is the 
thickness whose Tktb corresponds to the temperature 
of the experiment. (See Eq. (17) of Ref. ^3|.) By ad- 
justing d, they were able to obtain Q-AT data for both 
above and below the transition temperature. The pa- 
rameters that we obtained were d c = 5.4, z = 5.4 ± 0.4, 
and K g = 8.0 x 10~ 4 W/K. We assumed that K g was a 
constant over the parameter space that the Q-AT curves 
covered and estimated it based on the small d behavior of 
K. This placed some limits on this quantity but we were 
still able to vary it by a factor of two. We found that the 
scaling collapse was relatively insensitive to the value of 
Kg because most of the power was flowing through the 
film. As one can see, the scaling starts to break down for 
thicknesses (d ~ 9 layers) much larger than d c . 
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FIG. 4. The scaled Q-AT curves from superfluid 4 Helium 
films of (a) Ref. and (b) Ref. |54|. (Data set (b) has been 
shifted arbitrarily.) 

Curve b in Fig. || is from Ref. |54|. It contains more data 
corresponding to T > Tktb but doesn't cover as wida-a 
range of thicknesses overall as the data of Maps et alxJ 
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The collapse was found to be consistent with the results 
we have presented in this paper: d c = 5.2, z = 5.6 ± 0.5, 
and K g = 5.0 x 10~ 4 W/K. Like the JJA data of Ref. || 
the error bars here are large because the data extends 
only over a couple orders of magnitude. 

D. Universality 

One of the fundamental properties of critical behavior 
and scaling is universality: the idea that the same func- 
tion e±(x) in Eq. (|2l| ) and the same value of z describes 
all of the data independent of the system or material. For 
the T > Tktb branch, we found universality for nearly 
all of our scaled data. For the T < Tktb branch, on 
the other hand, the same function e±(x) described the 
data sets of the thinnest 2D samples but not that of the 
layered materialsc§C3 or the superfluid HeliumS We ex- 
plore these issues in this section. 

In Figure g the HTSC data from the YBCO juono- 
layer,o the LTSC data from the Hg/Xe thin filmj£3 and 
the YBCO JJA dataa are plotted together. Because the 
scaling functions are dimensionless, we multiply the x 
and y axes by non-universal constants which enter the 
scaling theory to account for the sample dependence. In 
order to test universality, all the data sets must be scaled 
with the same value of z. So we have adjusted each data 
set within the error bars of the parameters so that each 
has z — 5.6. As one can see, the collapse is very good, 
and bolsters the evidence that this data scales and. that 
z ~ 5.6 ± 0.3. Th£-.agreement of the Repaci et alE-i and 
the Herbert et aLcj data is particularly impressive since 
both data sets have extensive T < Tktb branches and 
because they come from two different systems: JJA's and 
SC's. 
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FIG. 5. n The JJA data of Fig. the YBCOjSC data 

of Fig. ^£3, and the Hg/Xe SC data of Fig. |cE3 plotted 
together. The scaled data sets have been shifted to show 
universality of the scaling function. 

In Figure ||, we plot all of the SC, JJA, and SF data 
that we have presented here in a single plot to test uni- 
versality. As in the previous figure, we adjusted each 
data set within the error bars so that z = 5.6. We also 
removed some of the individual scaled isotherms that did 
not scale well with the other scaled isotherms from the 
same sample. This, includes two of the isotherms from 
the Matsuda et alstB data (the second and third from the 
right in data set (a) of Fig. |^), parts of two of the data 
sets from data set (b) of Fig. || (second and third frcaa 
the right), and part of one data set from Harris et al$3 
(the rightmost set labeled (b) in Fig. |J). One can see 
clearly that the T > Tktb data scales well and strongly 
suggests universality for this temperature regime. The 
Garland and Leeo data is the weakest of these since it's 
lowest temperature isotherm does not lie on the other 
scaled data as we have pointed out in this figure. 
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FIG. 6. The data from all of the data sets plotted together. 
The T > Tktb data collapse very well suggesting universal- 
ity of the scaling function for that temperature regime. The 
T < Tktb data also scale well if three-data sets, belongingf-fcp 
layered materials (labeled AmmirataEpi and Vadlamannatil-3) 
and superfluid Helium (labeled Mapw), are neglected. This 
is discussed in the text. The Garland and Lee dataca scaled 
well except forrthe anomalous isotherm just above Tktb, la- 
beled Garland. H 

Universality was not found for all of the T < Tktb 
data. Besides the data plotted in Figure ||, there are 
only a few data sets that have significant T < Tktb 
branches: the HTSC BSCCO thin crystal of Ref. g|, the 
4 Helium data of Ref. |5^, and the YBCO multi-layer data 
of Ref. [l4| (We have labeled each of these data sets in 
Fig. ^.) The low temperature branches of these data sets 







fail to collapse with the other data. It is likely that the 
BSCCO and YBCO multi-layer data fail to collapse in 
this temperature range due to their layere d structure. 
We will discuss this point in Section IV A3 . As to why 
the T < Tktb 4 Helium scaled data do not collapse with 
the other data, a possible explanation is that the ap- 
proximate determination of £(T) was inaccurate. (Recall 
the lack of "ohmic" conductivity data for the superfluid 
measurements.) We recommend further studies of such 
data to more accurately determine £ (T) and to check the 
universality. In short, we cannot conclude whether the 
lack of universality of the scaling function for T < Tktb 
indicates a breakdown of the dynamic scaling results, or 
if it is due to the data or systems considered. We believe 
it is the latter. 



IV. DISCUSSION 



The main result of Section III is that z ~ 5.6 for super- 



conductors, superfluids, and Josephson junction arrays. 
The large value of z is nearly three times the expected 
value for these systems. Indeed there are many reports 
for superconductors and Josephson junction arrays that 
the TV exponent a jumps from 3 to 1 at the transition 
temperature, a result that is consistent with z — 2. In 
this section, we explore the discrepancy between the ap- 
proaches by discussing the validity of each. 



A. Dynamic Scaling 

Dynamic scaling is a powerful technique which has 
proved particularly useful in 

investigating the nature of the H-T phase diagram in 
the high-temperature superconductors. Yet, one must be 
careful with the results obtained with a scaling analysis 
because this techniq ue is not without its weaknesses. For 
example, in Section HI B , universality of the scaled data 
was examined and, while universality was found for the 
T > Tktb data and some of the T < Tktb data, it 
did not hold for all of the latter. In this Section, we will 
examine other aspects of the dynamic scaling analysis 
presented here and its validity. 



1. Vortex correlation length 

There are various aspects of assumptions that we made 
regarding the vortex correlation length that may detract 
from the dynamic scaling analysis. Here we make consid- 
erations on its symmetry, the value of the non-universal 
constant 6, and the expected range of the validity of the 
temperature dependence. 



Throughout Section [II it was assumed that the vor- 



tex correlation length was symmetric about the transi- 
tion temperature even though Ambegaokar et al.a predict 



that it should be asymmetric. [See Eq. (g).] Surpris- 
ingly, where there was data for both above and below 
the transition temperature, the scaling did not seem to 
suffer from this assumption. This is particularly so for 
the data on the superconductors. One notices that for 
the data of Refs. |2£| |34|, and ff3| the scaling works well 
for both branches. (See Figures ]l|and f§) 

When we allowed the correlation length to be asym- 
metric by allowing for different values of b below Tktb 
(&_), and above Tktb {b+), significantly better scaling 
could not be achieved, possibly indicating that £ is sym- 
metric, in accordance with the scaling of the numeri- 
cal results of Lee,-*|nd Teitel.c3 There were of course a 
few exceptions JfilrEJ but in neither case was prediction 
b + = 2nbJi verified. In fact, in these two cases, b- 
tended to be larger than b+. For the superfluid 4 He data 
of Ref. we found 6_ ~ 36 + . We do not view this 
result as conclusive because of the afore-mentioned prob- 
lems with determining £(T) for helium and we suggest 
further study of this topic. 

In our fitting of the resistance R(T) to determine 
£,±(T), we were sometimes able to fit the resistance to 
the Kosterlitz form over an extended temperature. range. 
A notable example is the YBCO mono-layer datacfJ where 
we fit the resistance from 25 K to 35 K. The upper limit of 
this fit is twice the Tktb which is remarkable since true 
KTB critical behavior is expected to be valid over a very 
narrow temperature range£ju3 This result could suggest 
that the critical region is larger than expected, or more 
likely that KT.B.-like behavior remains valid outside the 
critical regionEJ Indeed, Minnhagen and OlssonED indi- 
cate that Eq. (jlj) is a useful phenomenological form over 
wider temperature regimes as_Long as b is taken to be a 
phenomenological parametero We adopt this interpre- 
tation. 

In the literature, it is commonly-stated that b is mate- 
rial dependent but that b = 0(1). U In our results for su- 
perconductors, we found b to roughly O(10) but as small 
as 3.44 and 0.2, thereby varying by an order of magnitude 
or two. We also found that its value could vary within 
materials. For example, the value of b f or the 2 YBCO 
mono-layer systems examined in Section [II A varied by 
a factor of 4,0113 which is concerning. However, we do 
not think it is troubling that the b varied from mate- 
rial to material by an order of magnitude. The systems 
considered here are diverse from one another. For exam- 
ple, the electron density can vary significantly from the 
conventional superconductors to the HTSC's. Another 
interpretation is that the value of b is a phenomenologi- 
cal oru3-|ai}d not equivalent to its true asymptotic critical 
value .e3'E 



2. Universal jump: 6.6 to 1 

If one is to believe that z ~ 5.6 in these materials, then 
one expects to see a jump in the I-V exponent a of 6.6 to 
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1 in the J — > limit. Yet evidence of this is not observed 
in any of the samples. We believe this is because previous 
measurements have not gone to weak enough currents to 
observe this behavior. For example, based on the scaling 
curve of the YBCO mono-layer datao in Fig. [l]b where 
Tktb — 17.6K, one can see that the scaling curve is 
ohmic for x < 10 -5 . [x is the scaling variable.) This 
means that the 18 K isotherm would not become ohmic 
until the current I < 10~ n A. Clearly, voltage sensitiv- 
ity is far from detecting that crossover. [In fact, that 
isotherm would become first ohmic due to finite size ef- 
fects at a much higher current (Sec. IV filD -] It is a simi- 

taETw 



lar situation for the Herbert et al. datacd where the scal- 
ing curve is ohmic for x < 10~ 6 . The 74.722 K isotherm, 
the first isotherm above Tktb = 74.3 K, would not be- 
come ohmic until the current / < 10~ 8 A. This is several 
orders of magnitude smaller than the minimum current 
for that isotherm. 

The "jump" is evident in the scaled data only in the 
following way. For the T > Tktb data, the scaled data 
goes from ohmic behavior (a = 1) for small values of 
the scaling variable to a — 6.6 at large value of the scal- 
ing variable. In this way, the "jump" is only manifest 
as a crossover from the small x behavior to the large x 
behavior of t he T > Tktb branch of the scaled data. 

In Section IV C 2 , we will use a conventional approach 
to show that the behavior of a for the YBCO mono-layer 
datEO is consistent with the value of z that we found for 
that sample: z ~ 5.9. 



3. Three dimensional effects 

We have examined two samples which could be viewed 
as layered. (The first layered sample that we exam- 
ined was the YBCO/PBCO multi-layer system of Ref. || 
where the two unit-cell thick YBCO layers were believed 
to be electrically isolated. ' 
1000 A thick BSCCO crystal 



e second sample was the 
tB) In layered superconduc- 
tors, three dimensional (3D) behavior is expected in a 
small region of the current-temperature phase diagram 
near Tktb- One could then ask why such samples would 
scale in the same way as the much thinner samples and 
why 3D effects arc not manifest. 

To discuss the 3D effects, two new lengths are intro- 
duced. The first is related to the energy of the vortex 
pair. This length is the Josephson length A,/ = 7s, which 
incorporates the effects due to Josephson coupling be- 
tween the layers. (7 is the anisotropy ratio for layered 
superconductors and is equivalent to the ratio of the co- 
herence length in the ab planes to the coherence length in 
the c direction: 7 = £ a b/£c- s is the distance between lay- 
ers.) For separations less than Aj, the vortices interact 
with the 2D (logarithmic) potential. For larger separa- 
tions however, the potential becomes linear due to the 
Josephson coupling between the layersH (The Joseph- 
son coupling also introduces a term to the interaction 



energy for lengths less than Aj, but that term is very 
small compared to the 2D logarithmic one.) 

The second length is the c-axis vortex correlation 
length (to be distinguished from the c-axis super- 
conducting coherence length) and characterizes how far 
along the z direction the vortices are correlated. This 
length can also be viewed in terms of the length-scale de- 
pendent layer decoupling length £ 3D / 2D = 7^" (defined 
for T > T c ) since it determines the extent over which 
the 3D effects are important in the in-plane direction.^ 
Because of vortex screening, the Josephson interaction is 
screened out beyond lengths £^d/2Di making the interac- 
tion 2D at large separations for T > T c . £30/20 iiecomes 
small quickly above the transition temper ature.E3 

Because of these two competing lengths, 3D behav- 
ior is expected only over a small range of the TT phase 
space. On the temperature axis, this region is limited by 
s < £c(T) < D, where D is the thickness of the sample. 
On the I axis, the 3D region is limited to intermediate 



currents: Aj < r c < £ 



ZD 1 2D- 



Above ZWb, 



the renor- 



malized 7 (and hence Aj) grows quicklyo while £30 /2D: 
because of the temperature dependence of , decreases 
rapidly, further limiting the 3D behavior. This leaves 
only a small window in which 3D near Tktb effects can 
be observed. 

Returning our attention to the layered samples that 
we examined here, it is plausible that 3D effects are be- 
ing seen. After all, it is primarily the layered supercon- 
ductors of Refs. [44] and |43-that don't obey universality. 
For the multi-layer sampled, the T > Tktb scaled data 
fails to collapse with the other curves (see Fig. |^) for the 
isotherms nearest Tktb- That this sample could have 3D 
behavior is not in contradiction with the reports p£ those 
authors that their layers are electrically isolatedC-3 since 
magnetic coupling extends over larger distances. For the 
BSCCO sample of Ref. || whose T < Tktb data does 
not collapse with the others in Fig. ^, we have done the 
following calculation based on our above discussion to 
try to estimate where 3D effects should be seen. Since 
7 ~ 160 for BSCCO 2212, 3D effects should start to oc- 
cur for currents < 1.2 mA (where r c > Aj) and persist 
up to the current associated with the minimum of Lf s 
(~ 2.4mA) or the decoupling length, very near Tktb- 
Based on these rough estimates, which do not incorporate 
renormalization effects, it seems unlikely that 3D effects 
could be observed in these samples for T > Tktb- For 
the T < Tktb branch, it is more likely that the deviation 
from the universal curve is due to the thickness of that 
crystal since it is well known that the 3D region is rau< 
larger below the transition temperature than aboveES' 



B. Validity of "conventional" results 



In the previous Section ( IV A ) , we have addressed the 
validity of the dynamic scaling, whose results indicate 
that z ~ 5.6. Since this contradicts the evidence from the 
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conventional approach that z = 2, we now address those 
results. We do not claim that each paper is incorrect in 
their claims of z = 2 but we do hope to convince the 
reader that the question of the value of z is still an open 
one. 

The conventional results fall into roughly two classes, 
dc and ac. The dc measurements are the most com- 
mon and include the determination of a(T) from the I- 
V measurements, dc magneto-resistance measurements 
have also been used but less frequently. The ac measure- 
ments include the toKsiow oscillator 4 He measurements 
of Bishop and Reppyc2lcil and the ac inductance mea- 
surements of Fiory, Hebard, and Glaberson.tLa Clearly, 
we cannot address each paper that reports evidence for 
z = 2 and so we will discuss them in broad terms using 
particular examples where appropriate. 

1. I-V and R(T) dc Measurements 

Most of the papers that report evidence for a KTB 
transition or KTB behavior make their determinations 
based only on a(T) and R(T) measurements. (There 
are a rfew-. notable exceptions to this that we discuss 
below.t§c3) We suggest here that such an approach can- 
not be taken as evidence of KTB behavior and z = 2 but 
only as being consistent with z — 2 within the conven- 
tional approach. 

We begin with a brief description of this method. From 
the I-V data, a value of Tktb is determined from the 
condition, a(T = Tktb) = 3, (which, of course, as- 
sumes z — 2). It is then that this value of the transition 
temperatuEe-i is consistent with the R(T) data and the 
MinnhagerJHj form of the resistance: 

R(T) = Aexp[-2 v /fe(T c0 - T KTB )/{T - Tktb)]., (22) 

or the traditional form for R(T) [Eq. (jljj).] The mean 
field temperature T c q and constant A, in addition to b 
amount to three fitting parameters. (If Eq. ([to]) is used, 
then there are only two fitting parameters.) A further 
check that is sometimes used is to verify that a(T) de- 
creases quickly above Tktb- 

There are several reasons for why this approach can be 
a misleading check of z = 2. The first and most impor- 
tant is that determining where a(T) = 3 is difficult. It is 
well known that the predicted jump ^the-JJ/ exponent 
a(T) exists only in the I -> limitEaMHc3 and that, 
for temperatures above the transition temperature, the 
value of a(I,T) can vary quickly from a value of 1 to a 
value of z + 1 as a function of current. Hence, one must 
be sure that the value of reported a(T) will not dip to 
a lower value at currents whose voltages are below the 
voltage sensitivity. 

(Beyond the problem of probing the weak current limit, 
actually detecting the jump can be very difficult because, 
in the presence of finite size effects, small magnetic fields, 
and disorder, the jump gets rounded considerably, due 



to thepmami-pampeting length scales. Note that most 
papersQEllO'EJ do not report a jump in the I-V expo- 
nent.) 

The second problem with this approach is that, most 
of the time, a misleading criterion is used to determine 
a(T). Because it is seldom (if ever) the situation that one 
knows that the value of a(T) determined for a particular 
isotherm represents the weak current limit, one should 
examine a(T) for a given length scale (or cpmmon cur- 
rent), as originally suggested by Kadin et alE5 Here, one 
looks for a change in the behavior of this quantity near 
the transition temperature. This approach is usually not 
followed rigorously however. Instead of using a common 
current, investigators will determine a(T) for a common 
voltage range. This has the effect of biasing the results 
because, as one looks at higher temperature isotherms at 
common voltages, one is looking at lower currents. This 
means that a(T) will have decreased not only due to tem- 
perature but also due to longer length scales. Therefore, 
any report rapid decreases of ot(T) may be an artifact of 
examining the isotherms at a constant voltage. (There 
are also cases where not even a common voltage range is 
used; rather, the parts of the isotherms that deliver the 
desired rapid change in a(T) are studied.) 

Thirdly, once a value of Tktb is obtained for where 
a = 3, the number of fitting parameters in Eq. (|2^ ) 
used to verify the self consistency is large (three) and 
the temperature range over which one fits the R(T) data 
is limited. So we contend that the conditions for check- 
ing z = 2 in this approach are not stringent enough to 
be considered as evidence. 

Finally, we note the inherent problem with this ap- 
proach mentioned above that z is not allowed to vary in 
order to optimize the analysis. 

As an example, in Ref. |49|a jump in the I-V character- 
istics is reported for the YBCO/Ag weak link J J A. (See 
Fig. 2 inset in Ref. ||a or Figs. 5-7 and 5-8 in Ref. ||b.) 
The sharp jump in Ref. ^ was obtained by fitting the 
ohmic part of the curves at a temperature just above 
where the I-V has a = 3. In this case, no fit of R(T) was 
done to check the value of Tktb for consistency. Also, 
one cannot rule out that, if those authors had another 
decade or two of voltage sensitivity, isotherms at temper- 
atures lower than their Tktb would also become ohmic 
at lower currents. Further, as is clear from the scaling 
analysis shown in Figs. [|a and [|, this data is consistent 
with z = 5.6. As a second example of this, we look at the 
YBCO mono-layer data of Ref. |45|. By inspecting Fig. 4 
of the reference, it is clear that some of the isotherms be- 
low their reported Tktb could manifest ohmic behavior 
if more decades of voltage sensitivity were possible. An 
inspection of the I-V data of Ref. |44| (Fig. 2) yields a 
similar conclusion. It is unlikely that the reported value 
of a — 3 at T = Tktb = 40.1 is the asymptotic value of 
that quantity at low currents. [The voltage sensitivity for 
that data also occurs at a larger value (~ 10~ 7 V) than 
that of other measurements (~ 10 -9 V).] 

To summarize, I-V along with R(T) measurements do 
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provide self-consistent results for a(T = Tktb) = 3 (and 
thus z = 2), but cannot be taken as proof that z = 2. 
This is because the flexibility in determining; Tkt r and 
the three fitting parameters used in Eq. (22) to fit the 
smooth, monotonic R(T) data do not pose tight enough 
constraints. 



2. Kinetic inductance on SC's, magneto-resistance, and 
Helium torsion experiments 



We now turn to other measurements used to study 
dynamics in the SC, SF, or JJA systems. 

As mentioned above, there are a few notable pa- 
pers that went beyond measuring only I-V and R(T) 
curves. The most comprehensive study was done by 
Fiory, Hebard and GlabersonEa who looked at kinetic in- 
ductance and magnetoristance in addition to the usual 
I-V and R(T) measurements on In/InO systems. The 
thoroughness of their approach is commendable. While 
their results seem to be consistent with z = 2, they are 
not without their inconsistencies and cannot necessarily 
be taken as definitive evidence for z = 2. 

Fiory et alr3 were able to determine the vortex in- 
teraction strength q 2 in two ways for a rectangular strip 
sample of In/InO. In their ac impedance measurements 
used to determine the kinetic inductance, the frequencies 
are small (160 Hz) and so they are probing only stat- 
ics. In that case the kinetic inductance does not depend 
upon the density of free vortices rif but only the super- 
fluid density. In that case, q 2 is directly determined. In 
their I-V measurements, a dynamic value of q 2 is deter- 
mined. They find that the values of q 2 from the two 
measurements agree over a limited temperature range 
(~ 30mK), which would be consistent with a value of 
z = 2. (See Fig. 9 of Ref. However, this result 

should be viewed cautiously since the two measurements 
disagree for most of the region T < Tktb- Further, the 
four-probe contact method that they used to measure the 
data for the comparison of the values of q 2 is less accurate 
and less sensitiveEa than the two-coil contactless method 
that they used on circular sample from the same film. 
Moreover, their measurements of samples from the same 
film revealed variations in T c that suggest that sample 
inhomogeneities may be significant. 

To further-test the value of Tktb for the In/InO films, 
Fiory et aZ.lla measured the magneto-resistance R(H). 
According to the theory of Minnhagenj£3 R(H) should 
be linear in H at Tktb but sub-linear above it and 
faster than linear below it. Their data at T = 1.782K 
do show a region which is nearly linear (R cx iJ 107 ) over 
roughly two decades. This temperature for crossover is 
roughly consistent with their Tktb — 1.782K determined 
from a(T = Tktb) = 3 criteria. As mentioned above, 
the roughly linear area is over only two decades and the 
samples do have a degree of inhomogeneity to them. Fur- 
ther, it is4ikcly that surface barriers should be taken into 
account E2 



Garland and Lee@ have also used magneto-resistance 
data in addition to the I-V and R(T) measurements 
on In/InO films. Based on their Minnhagen criteria for 
R{H) they determine Tktb = 3.123K. At this tempera- 
ture, they find that a(T) drops from a value of roughly 
four to one. One will also notice that their I-V at Tktb is 
not a pure power law as required by the dynamic scaling 
and that their T < Tktb isotherms all have a positive 
curvature. They attribute this behavior to field-induced 
vortices. It is clear from Fig. 7 in their paper that the 
field plays a role for fields at least as low as 5 mG. Yet, in 
keeping with the discussion on finite size effects in Sec- 
it is our view that the crossover to ohmic 
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tion 

I-V in their Fig. 5 should occur at a roughly common 
value of the current. This is because the magnetic length 
scale should not change significantly over the 100 mK 
that their I-V isotherms cover below their claimed Tktb- 
This is not observed in that data. Choosing Tktb at a 
lower temperature seems to be a better explanation, es- 
pecially when one considers the good collapse of their 
data in the dynamic scaling analysis as seen in Figure 



4 Helium torsion experimenta23'Ell were among the first 
evidence for the Kosterlitz-Thouless-Berezinskii transi- 
tion. Nevertheless, while we agree that these measure- 
ment are indicative of KTB behavior (i.e., vortex pair 
unbinding), we do not believe that they indicate z — 2. 
In the case of 4 Helium torsion experimentsEJ, Bishop and 
Reppy measured the period shift and Q value of an An- 
dronikashvili cell. The former is proportional to the real 
part of the dielectric constant e and the latter to the 
imaginary part of that quantity. It is the latter that 
has the predominant dependence on the free vortex den- 
sity n f . [See their Eq. (A2) or Eqs. (3.17) and (3.18) 
in Ref. |2|.] They implicitly assume z = 2 in writing 
Im(e) oc nf. We point out that their method cannot 
distinguish a value of z but only the product bz. A value 
of z other than 2 would still make the values of their 
fitting parameters reasonable. The noise spectrum mea- 
surements of Shaw et a/J13 are another good example of 
a measurement that is not lable to determine the value of 
z but only the product bzEJ Hence these measurements 
cannot be taken as evidence for z = 2. 



C. Comparison of Conventional Approach with 
Dynamic scaling 

To further examine the validity of the finding z ~ 
5.6 ± 0.3 obtained from the Fisher-Fisher-Huse scaling, 
we will examine a data set using the conventional ap- 
proach with an arbitrary value of z. We have chosen 



the YBCO mono-layer data of Ref. 34 for several rea- 
sons. First, their data covers by far the largest current, 
voltage, and temperature range of any data set in the 
literature. Secondly, those authors have pointed out that 
there data does not satisfy the "conventional" criteria for 
KTB behavior: the I-V exponent does not vary rapidly 
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from 3 to 1 near the transition temperature. Finally, they 
suggested finite size effects to explain the lack of KTB 
behavior and performed a largely qualitative analysis. 

In this section, we will use a quantitative analysis to 
show that finite size effects cannot account for the ob- 



served behavior in that data (Sec. [TV C l| .) We will then 
perform a conventional analysis of the I -V exponent to 
look for evidence of z ~ 5.9 (Sec. |IVC2 .) 



1. Finite size effects 

Repaci et a/S have been suggested that their YBCO 
mono-layer data that has been shown here to scale so well 
with z = 5.9 can be explained in a "conventional" man- 
ner by finite size effects, following similar suggestions re- 
garding low current ohmic behavior of others .EjoS We 
investigate that possibility in this subsection, making an 
explicit comparison between the dynamic scaling analysis 
and the conventional approach. As the reader will see, 
we find no evidence of finite size effects in any of the data 
that we examine. 

This subsection is organized as follows. After a dis- 
cussion of the principle differences between the two sce- 
narios, we look for evidence of finite size effects in the 
I-V data from the YBCO mono-layer dataE.3 first using 
the conventional picture and then in the dynamic scaling 
analysis. 

Theory 

The primary difference between the finite size expla- 
nation and the dynamic scaling explanation of that data 
is the following: the former ascribes finite-size induced, 
low current ohmic behavior to I-V isotherms with tem- 
peratures T < Tktb', in the dynamic scaling explana- 
tion the low current ohmic behavior is associated with 
T > Tktb isotherms. The issue of course is the place- 
ment of the transition temperature. Because the scaling 
analysis indicates that u{Tktb) — 6.8 and the conven- 
tional analysis assumes C((Tktb) = 3, Tktb will be at a 
lower temperature in the scaling analysis scenario. 

We first examine the nature of the ohmic to non-ohmic 
crossover in the finite-size-effect scenario. Recall that 
with finite size effects, there is always a density of free 
vortices. [See Eq. (||).] This means that at small currents 
where one is probing large length scales, the free vortices 
will dominate the resistance and the I-V characteristics 
will be ohmic. As one probes smaller length scales by in- 
creasing the current (r c <C Lfs), current-induced vortex 
unbinding will start to dominate the resistance and the 
I-V curves will become non-ohmicS Because either the 
system size or the 2D penetration depth can be on the 
order of 1-10 microns, one would expect finite size effects 
to have an influence. 

We will now use the conventional picture to make the 
above discussion more quantitative for T < Tktb- We 
will show that the I-V curves are ohmic when r c ^> Lf s 
and non-ohmic when r c -C Lj s . With finite size effects, 



the energy of a vortex pair is generalized from Eq. (g) 
to a potential that is logarithmic^r R <c Lf s and ap- 
proaching a constant for R ^> L / S -E3 To approximate this 
behavior, we use 

E(R) = q 2 {ln(R/Zo) - HR/Lfs + 1)] - JR/Zo + 2E C 

(23) 

where q 2 = itn 2D Ti 2 /2m and J = Trhl^d/eA is a cur- 
rent with dimensions of energy. It is the second term on 
the right hand side of the equation that causes the zero- 
current (J = 0) pair energy to approach a constant as 
1/R for R> Lf s . With this term, E(R) no longer peaks 
at r c . Rather, it peaks at 



Js 



Lfs_ 
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LfsCoQ 2 /J 



(24) 



There are two limits to this equation. The large cur- 
rent limit J 3> q 2 £,o/Lf s can be rewritten as r c <C Lf s 
while the small current limit can be expressed, r c 3> 
Lf s . In the large current limit where one is probing 
small lengths, r[ s approaches r c while at small currents 



( J < q 2 £o/L fs 



-J Lf s ^oq 2 j J . Remember that it 



is the value of E(R) at r[ s (i.e., the height of the barrier) 
that determines the density of free vortices. In one limit 
(r c « L fs ), 



E(rl s )=q 2 [Hq 2 /J)-l}. 
On the other hand, 

E(r[ s )^q 2 HL fs /Co) 



(25) 



(26) 



for r c » Lf s . So, for T < Tktb, the I-V's are ohmic at 
small currents since E(r[ s ) does not depend on current. 
At large currents, E(r[ s ) depends upon the current in the 
traditional way, and one finds the usual I-V relationship: 
V cx I a . The ohmic to non-ohmic crossover occurs when 
/ — Ak B T K TBWe/iTfiLfS (i.e., r c ~ Lf„.) 

The T > Tktb effect of finite size on the transport 
behavior is a little more subtle than the low temperature 
effect because another competing length scale is involved. 
Even in the absence of finite size effects for this temper- 
ature range, the I-V curves cross over from thermally 
dominated ohmic behavior at small / to current-induced 
non-ohmic behavior at large /. The current at which this 
crossover occurs depends upon the size of the largest vor- 
tex pairs, ^+(T),E3 a quantity which is strongly tempera- 
ture dependent. With finite size effects, since all currents 
that probe lengths greater than L f s are ohmic, L f s com- 
petes with £+(T) yielding the following conditions: when 
r c 3> min[L/ s , £+] there is ohmic behavior and non-ohmic 
when r c <C min[X/ s , £_)_]. An important term in our dis- 
cussion is "premature-ohmic" behavior, which denotes 
ohmic behavior at currents for which there would not be 
ohmic behavior in the absence of finite size effects. It 
occurs when Lf s < (and when the temperature is 

sufficiently close to Tktb-) 
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It is the temperature dependence of these conditions 
(as well as the magnitude of the crossover current) that 
mark the signature of finite size effects for T > Tktb- 
For T > Tkt b , £+ (T) is large and exceeds Lj s so that the 
value of the current at which the isotherms cross from one 
behavior to another will depend upon Lf s and not the 
size of the largest pairs. In this case, the crossover cur- 
rent will depend only weakly on temperature. As £+(T) 
becomes smaller than L f s at temperatures further above 
the transition, it is that sets sets the current scale 

for the ohmic to non-ohmic crossover. So the crossover 
current becomes strongly temperature dependent. 
Conventional check of finite size effects 

As mentioned above possible candidate for observing 
finite size effects is the mono-layer YBCO data of Ref. |34| 
because the X2D — 40/xm. To determine the value of 
the current at which "premature-ohmic" behavior should 
occur for this data, we use r c = Lf s = AO/im. to solve for 
the current. Using W = 200/im, we find I cr0 ssover — 
2.4/iA, a value much smaller than the observed crossover 
current: ~ 100/xA for T ~ 28K. This would seem to 
indicate that A213 is not responsible for the ohmic to non- 
ohmic crossover in these materials. A renormalization 
group studju confirms that this remains the case after 
renormalization effects of Lj s are accounted for. This is 
because the condition for observing the finite size effects 
is r c > Lf s . Under renormalization, Iu~ does shrink 
quickly but r c shrinks even more quickly.119 

The length scale that corresponds_to the approximate 
crossover current for the Repaci et alEj data for tempera- 
tures around 28K is l//m. Because neither a field-induced 
vortex length or a pinning (disorder) length could corre- 
spond to this value, it seems unlikely that the behavior 
observed in Ref. [54] is due to finite size effects. 

Not only is the magnitude of the crossover current in- 
consistent with finite size effect but so also is the temper- 
ature dependence of this quantity. If it were L f s and not 
which determined the ohmic to non-ohmic crossover, 
the data would not collapse so well in the dynamic analy- 
sis since the temperature dependence of these two quan- 
tities are so different. 

Another signature for finite size effects in the con- 
ventional picture is that the resistance will have 
an arrhenious temperature dependence: R(T) = 
Ajy exp{B/KBT} where B is related to the energy of a 
free vortex and An is a constant. [See Eq. (B).] We have 
fit this formula to the data, as shown in Fig. || (dotted 
line), but do not get a satisfactory result. .— ■ 

Finally, a finite size analysis of Repaci et al£n data 
would indicate a Tktb in excess of-,30K following the 
subsequent work of Herbert et alS3 This is contrary 
to the mutual inductance data of Gasparov et aZ.O on 
YBCO mono-layer films which probes the temperature 
at which the largest pairs unbind. (Remember, the size 
of the largest pairs is expected to decrease as £+ (T) above 
Tktb-) They find that vortex pairs of size ~ 0.018/im 
unbind at a temperature of roughly 47. OK, vortex pairs of 
size ~ 0.78/im unbind at a temperature of roughly 27. 8K, 



and vortex pairs of size ~ 1.58/zm unbind at a tempera- 
ture of roughly 25. 5K. This trend is consistent with the 
value of Tktb that we find (17. 6K) for the YBCO mono- 
layer data of Ref. |34| using a dynamic scaling analysis. 
Dynamic scaling check of finite size effects 

The signature of finite size effects in the dynamic scal- 
ing depends upon the temperature. For T < Tktb, 
the scaled data should peal off the scaling curve to 
go ohmic (slope 1 in Figs. |l[|^) for currents less than 
4kBTKTBWe/irh~L f s . For temperatures above the tran- 
sition temperature, one would observe finite size effects 
only if L f s were shorter than the size of the largest vor- 
tex pairs. And in that case, the scaled data would break 
from the scaling curve to become "prematurely-ohmic" 
at I = AkBTKTBWe/irfiLfs. This is not observed in 
any of the scaled data in this paper. Such behavior was 
observed however for the BSCCO data of Ref. E. (See 
Fig. 2 there.) In that case however, it seems more likely 
that this behavior is not due to finite size effects but 
to voltage sensitivity. The crossover to the premature- 
ohmic behavior occurs more rapidly than one would ex- 
pect for finite size effects and also occurs at roughly the 
same voltage, which is near the voltage sensitivity limit. 
Finite size effect discussion 

As mentioned above, there is not a "true" thermody- 
namic phase transition in superconductors because of the 
finite penetration depth. One could then ask why there 
is a critical isotherm at all. For example, in the Repaci 
et alE3 data, the 17K isotherm is straight over nearly 9 
decades of voltage. The answer to this is that the cor- 
relation length at finite current does not become infinite 
(and is not longer than Lf s ) even though it is very close 
to T c . This is apparent from Eq. (||) where it is seen that 
the correlation length decreases as the reciprocal of cur- 
rent. It is only when £(/, T) ~ Lf s that one would begin 
to see deviations in the critical isotherm. 

To summarize this section, we find no evidence of fi- 
nite size effects in any of the data that we examine. In 
principle however, finite size effects are inherent to su- 
perconductors and will manifest themselves if the prob- 
ing current is small enough. It is our opinion that none 
of the data sets that we examined went to currents small 
enough to detect finite size effects. 

2. Conventional Approach 

We now examine the Repaci et al. data with a con- 
ventional approach. In Figure [?], we plot the I-V ex- 
ponent (V oc 7 Q ( T )) at a fixed current (I = 0.7mA) as 
determined from the <f[log V]/d[\ogI] data in Figure 3 of 
Ref. HI). (The error bars were determined from that figure 
and from fits to the TV curves, and are only shown for 
the near-linear region. The a(T < 17K) data also came 
from power-law fits to the digitized I-V data from Fig. 2 
of that reference.) As one can see, there are no features at 
a ~ 3 that would suggest a phase transition, as originally 
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pointed out by those authors. An interpretation of this 
data is difficult. A possible feature is a crossover from 
near-linear behavior of a(T) to nonlinear behavior near 
the value of Tktb (17.6K), obtained from the dynamic 
scaling behavior. Further, the value of a observed at this 
temperature produces an estimate for z (~ 7) that is sim- 
ilar to that obtained from the scaling procedure. One can 
see that the T > Tktb behavior is concave up, contrary 
to the analytical work of Ref. ^8], but more consistent 
with the simulational work of Refs. \TA 
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FIG. 7. a(T) data at a fixed current (I = 0.7mA) for an 
YBCO mono-layer taken from the d[log V]/d[log I] data in 
Figure 3 of Ref. Also plotted are a linear fit to the data 
less than the transition temperature and the lines a(T) — 7.1 
and a(T) = 1.0 used to determine the parameters T c q, z, and 
to- [INSET: a(T) data at a fixed current (J ~ 4.0mA) for the 
BSCCO film of Ref. ||] 

The values of Tktb and T c o determined from the stan- 
dard analysis of a(T) with arbitrary z are less definitive 
because one is not assuming a value of a(T^Ts)- To 
be consistent with the value of Tktb determined from 
the scaling analysis, we have chosen Tktb = 17. 6K. The 
subsequent value of T c q determined from a linear fit to 
the a(T < 17.6K) is T c0 = 24.8K. (See Figure @.) z was 
found to be 6.1 ± 0.2, consistent with the dynamic scal- 
ing value. The renormalized dielectric constant also has 
a reasonable value: e c = 1.59. 

One possible explanation for the absence of a clear 
jump in the exponent a(T) in Figure is the relatively 
short length scale it represents. In order to obtain well- 
defined and stable a(T) values (at the same reference 
current for all of the isotherms), it was necessary to 
choose a relatively large value for the reference current. 
This current value corresponds to the length scale 1400A 
(<C Lf s ), which is distant from the desired I — ► limit. 
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FIG. 8. The resistance data of Ref. and fits to Eq. (|T|) 
(dashed line) and Eq. (^) (dotted line). The latter equation, 
based on assumptions of finite size effects, does not adequately 
describe the data. 

The final step in the conventional approach is to ex- 
amine R(T). Eq. Ji"s| ) was used in place of Eq. (22) to 
fit the data since many of the isotherms are obtained in 
the regime T > T c0 = 24. 8K, where Eq. @ is not valid. 
The dashed line in Figure |§| is the fit using the parameters 
that optimized the scaling for Fig. [lb: Tktb — 17.6 K, 
b = 7.79, and z = 5.9±0.3. Clearly, the data is more con- 
sistent with KTB critical behavior than finite size effects 
discussed in the previous section. 

We repeat this analysis for the BSCCO sample (1000A 
thick) of Ref. [l3| The a(T) determined for a constant 
current is shown in the inset of Fig. |?]. As one can see, 
this data is more noisy and covers a much smaller temper- 
ature range than that in Fig. 0, thereby precluding a com- 
plete conventional analysis. So, while the value of a at 
Tktb = 78.87K is consistent with that determined from 
the dynamic scaling analysis (6.6) and a seems to change 
behavior at that temperature, one could not claim these 
observations as evidence for u(Tktb) = 6.6. However, 
an important observation can be made by comparing this 
inset to the inset of Fig. 1 of Ref. [5^, which shows a(T) 
determined from the same data, but for a constant volt- 
age. As one can see, a(T) decreases much more rapidly 
for a constant voltage th an for a constant current, rein- 
forcing our claim in Sec. [VB 1 that a rapid decrease in 
a(T) could be an artifact of using the constant voltage 
a(T) data. 



D. Theoretical Considerations 

The primary degrees of freedom associated with the 
KTB phase transition are vortices. The dynamic behav- 
ior should therefore be dissipative. Specifically, it has 
been arguedli^ that superconducting dynamics in zero 
field may be purely relaxational (model Ac3) for any di- 
mension, with a diffusive exponent, z < 2. This inter- 
pretation is consistent with the conventional treatment 
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of KTB dynamics. However, the present scaling analysis 
of I-V and Q-AT data from SC's, JJA's and SF's indi- 
cates that z ~ 5.6, a result consistent with sub-diffusive 
dynamics. Here we mention some possible explanations 
for this large value of the dynamic exponent. 

Pinning is known to play a crucial role in the large 
values of z £s> 2\ observed in vortex glass phenomena in 
experimentaI3~L3 and Monte Carlo simulationsBuO in 
high temperature superconductors. However, we do not 
believe that pinning can explain the surprising values of 
z obtained in the present 2D analysis, at zero field. The 
reason is that our result, z ~ 5.6 is obtained from very 
distinct systems: superconductors, JJA's, and superffu- 
ids. For superfluids, in particular, pinning effects should 
be negligible. A pinning explanation therefore appears 
inconsistent with the universal nature of z. _ 

Collective excitations, such as vortex density waves,L3 
may mediate the observed dynamic behavior. This is 
more likely to be true if the vortices cannot exit the sam- 
ple easily, because of surface barriers. We can determine 
the dynamical critical exponent for this behavior as fol- 
lows. Based upon the Coulomb gas analogy, the vortex 
plasma frequency is given by ui cx y/n, where n ~ 
is the vortex density. Using t oc we find z — 1.5 in 
three dimensions and z — 1.0 in two dimensions. Vortex 
plasmons therefore cannot explain the large values of z. 

Another possibility is that the suppositions leading to 
Eq. ([l7| ) are incorrect. To arrive at that equation, only 
two forces are included, a viscous force and the Lorentz 
force. Perhaps, with -the inclusion of other forces (e.g. 
surface barrier forces,E3) an explanation for z — 5.6 can 
be found. 

We believe the most likely explanation for large values 
of z lies in copulated vortex motion, described as. "part- 
ner transfer, "E3 or "collaborative dissociation" c3 These 
suggest mechanisms whereby bound vortex pairs do not 
simply dissociate into free vortices. Instead, the pro- 
cess is mediated by neighboring vortex pairs, in terms of 
consecutive recombination-dissociation events. Further 
work is required to confirm this model. 



V. SUMMARY 

As stated in the introductory paragraph, the dynamic 
scaling approach presented here is different than most 
previous studies of dynamics in 2D SC's, JJA's, and SF's 
in that this approach allows one to vary z to optimize the 
analysis. In the "conventional" approach, the value of z 
is implicitly taken to be two. By using the dynamic scal- 
ing analysis and allowing the value of z to vary, we have 
presented evidence which suggests non-diffusive behav- 
ior. 

Via a dynamic scaling analysis of transport data from 
SC's, JJA's, and SF's, we find z ~ 5.6 ± 0.3, contrary 
to the value assumed but not tested in previous reports. 
This analysis seems convincing in that the collapse is 



excellent in many data sets and the value of z is robust 
from system to system and material to material. 

The results of the dynamic scaling analysis also go 
against the many studies consistent with z = 2. We 
have included in this work a discussion of those "conven- 
tional" approaches to studying the dynamics of the KTB 
transition. Like the dynamic scaling analysis, we find 
that this approach also has its drawbacks. Perhaps the 
most important is that the experiments do not yet have 
the sensitivity to actually observe the predicted jump in 
the I-V's in the I — ► limit which-,we estimated to be 
at 10~ n A for a particular sample£-3 Another drawback 
is that these approaches do not vary z to optimize the 
fits. Furthermore, the most common method of "verify- 
ing" KTB behavior [which is to obtain Tktb from the 
condition oi(Tktb) = 3 and then to do a three parameter 
fit to Eq. (|22|),1 does not pose constraints tight enough 
to prove z = 2. Our analysis of the evidence for z = 2 
raises the following questions: 

• Why has no scaling of zero-field I-V data with z = 
2 been realized; 

• Why do critical isotherms have a much larger value 
of a than the value consistent with z = 2; 

• If finite size effects are present, why does the ohmic 
to non-ohmic crossover not coincide with r c = Lf s . 

We also compared directly the conventional approach 
and the dynamic scaling approach for data from a par- 
ticular sample: the YBCO mono-layer data of Ref. |34|. 
In that reference, those authors found that the conven- 
tional approach can not explain their data. In Section 
IV C 1 , we found that an incorporation of the finite size 



effects into the conventional approach is also not consis- 
tent with their data. Further in Section IV C 2 . we saw 



that a conventional analysis of the I-V exponent was con- 
sistent with z ~ 5.6. A dynamic scaling analysis of their 
data however resulted in a beautiful collapse, as shown 
in Fig. 0b. 

The primary purpose of this paper is to convey that 
the question of the value of z in these systems is still an 
open one, despite the conventional wisdom that 2 = 2. 
We believe that more study is needed. In particular, more 
data on all systems, especially JJA and SF is needed and 
over wider temperatures and current regions (or Q r. 
gions for SF's). The impressive data of Repaci et al. 
sets a good standard. Not only should dynamic scaling 
analysis be tried on this data but so too should compre- 
hensive "conventional" studies like those of Fiory, Hebard 
and Glaberson.E3 (By "comprehensive" , we mean going 
beyond just the usual o.{Tktb) = 3 and R(T) measure- 
ments.) Of special importance would be a measurement 
of q 2 using static kinetic inductance data (in the appro- 
priate frequency range) and the dynamic I-V exponent 
a. If the conventional theories are valid after general- 
ization to a general z, then a — 1 = [z/2]q 2 . (Fiory et 
al. did do such a measurement but found agreement only 
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over 30mK.) Allowing the value of z in the conventional 
theory would also be a useful exercise. 

The authors gratefully acknowledge conversations with 
S. M. Girvin, L. I. Glazman, B. I. Halperin, A. Hebard, 
J. M. Kosterlitz, C. J. Lobb, P. Minnhagen, P. Muzikar, 
and S. Teitel. This work was supported by the Midwest 
Superconductivity Consortium through D.O.E. Contract 
No. DE-FG02-90ER45427 and by NSF Grant No. DMR 
95-01272. Acknowlegement (SWP) is made to the donors 
of The Petroleum Research Fund, administered by the 
ACS, for support of this research. 



* Electronic address: pierson@wpi.edu. 
' Present 

address: Applied Superconductivity Center, University of 
Wisconsin-Madison, 1500 Engineering Drive, Madison, WI 
53706. Electronic address: mgfriesen@facstaff.wisc.edu. 

* Present address: 

1 V. Ambegaokar, B. I. Halperin, D. R. Nelson and E. D. Sig- 
gia, Phys. Rev. Lett. 40, 783 (1978). 

2 V. Ambegaokar, B. I. Halperin, D. R. Nelson and E. D. Sig- 
gia, Phys. Rev. B 21, 1806 (1980). 

3 See, e.g., M. Takahashi, J. Phys. Soc. Japan 48, 746 (1980). 

4 See e.g., F. G. Mertens, A. R. Bishop, G. M. Wysin, and 
C. Kawabata, Phys. Rev. Lett. 59, 117 (1987). 

5 See, e.g., A. Jonsson and P. Minnhagen, Phys. Rev. Lett. 
55, 9035 (1997); and references therein. 

6 B. J. Kim, P. Min nhagen, and P. Olss on, unpublished, Los 
Alamos Archives: 



cond-mat/9806231 



7 J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 
(1973). 

8 J. M. Kosterlitz, J. Phys. C 7, 1046 (1974). 

9 V. L. Berezinskii, Sov. Phys. JETP 32, 493 (1971). 

10 D. S. Fisher, M. P. A. Fisher, and D. A. Huse, Phys. Rev. 
B 43, 130 (1991). 

11 S. A. Wolf, D. U. Gubser, and Y. Imry, Phys. Rev. Lett. 
42, 324 (1979). 

12 D. C. Harris, S. T. Herbert, D. Stroud, and J. C. Garland, 
Phys. Rev. Lett. 67, 3606 (1991). 

13 L. Miu, P. Wagner, U. Frey, A. Hadish, D. Miu, and 
H. Adrian, Phys. Rev. B 52, 4553 (1995). 

14 R. Brada, H. Chayet, and W. I. Glaberson, Phys. Rev. B 
48, 12874 (1993). 

15 R. Theron, J.-B. Simond, Ch. Leemann, H. Beck, and 
P. Marinoli, Phys. Rev. Lett. 71, 1246 (1993). 

16 D. Bormann, Phys. Rev. Lett. 78, 4324 (1997). 

17 P. Minnhagen et al, Phys. Rev. Lett. 74, 3672 (1995); see 
also Ref. hq, and references therein. 

18 A. T. Fiory, A. F. Hebard, and W. I. Glaberson, Phys. Rev. 
B 28, 5075 (1983). 

19 T. K. Shaw, M. J. Ferrari, L. L. Sohn, D.-H. Lee, M. Tin- 
kham, and J. Clarke, Phys. Rev. Lett. 76, 2551 (1996). 

20 D. J. Bishop and J. D. Reppy, Phys. Rev. Lett. 40, 1727 
(1978). 

21 D. J. Bishop and J. D. Reppy, Phys. Rev. B 22, 5171 



(1980). 

22 B. I. Halperin, in Physics of Low-Dimensional Sys- 
tems, Proceedings of the Kyoto Summer Institute, 1979, 
edited by Y. Nagaoka and S. Hikami (Publication Office, 
Prog. Theor. Phys., Kyoto, 1979) p. 53. 

23 J. E. Mooij, in Percolation, Localization and Superconduc- 
tivity . Proceedings of a NATO Advanced Study Institute, 
edited by A. M. Goldman and S. A. Wolf, (Plenum; New 
York, NY, USA 1984) pp. 325-70. 

24 P. M. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987). 

2o A. Hebard, in Strongly Correlated Electronic Materials: 
The Los Alamos Symposium, 1993, Edited by K. S. Be- 
dell, (Addison Wesley, New York, 1994), p. 251. 

26 M. V. Simkin and J. M. Kosterlitz, Phys. Rev. B 55, 11646 
(1997). 

27 J. Pearl, Appl. Phys. Lett. 5, 65 (1964). 

28 A. M. Kadin, K. Epstein, and A. M. Goldman, Phys. 
Rev. B 27, 6991 (1983). 

29 B. I. Halperin and D. R. Nelson, J. Low Temp. Phys. 36, 
599 (1979). 

30 V. Ambegaokar and S. Teitel, Phys. Rev. B 19, 1667 
(1979). 

31 M. R. Beasley, J. E. Mooij, and T. P. Orlando, Phys. Rev. 
Lett. 42, 1165 (1979). 

32 S. T. Herbert, Y. Jun, R. S. Newrock, C. J. Lobb, K. Ravin- 
dran, H.-K. Shin, D. B. Mast, and S. Elhamri, Phys. Rev. 
B 57, 1154 (1998). 

33 See C. J. Lobb, D. W. Abraham, and M. Tinkham, Phys. 
Rev. B 27, 150 (1983), and references therein. 

34 J. M. Repaci, C. Kwon, Q. Li, X. Jiang, T. Venkatessan, 
R. E. Glover, C. J. Lobb, and R. S. Newrock, Phys. Rev. 
B 54, R9674 (1996). 

35 D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39, 
1201 (1977). 

36 S. W. Pierson, Phil. Mag. B 76, 715 (1997). 

37 B. A. Huberman, R. J. Myerson, and S. Doniach, Phys. 
Rev. Lett. 40, 780 (1978). 

38 J. Bardeen and M. J. Stephen Phys. Rev. 140, A1197 
(1965). 

39 See, N. Goldenfeld, Lectures on Phase Transitions and the 
Renormalization Group (Addison- Wesley, Reading, MA, 
1992), p. 213. 

40 P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49, 
435 (1977). 

41 J. Lidmar, M. Wallin, C. Wengel, S. M. Girvin, and 
A. P. Young, Phys. Rev. B 58, 2827 (1998). 

42 J.-R. Lee and S. Teitel, Phys. Rev. B 50, 3149 (1994). 

43 S. M. Ammirata, M. Friesen, S. W. Pierson, L. A. Gorham, 
J. C. Hunnicutt, M. L. Trawick, and C. D. Keener, Physica 
C (in press). 

44 S. Vadlamannati Q. Li, T. Venkatesan, W. L. McLean, and 
P. Lindenfeld, Phys. Rev. B 44, 7094 (1991). 

45 Y. Matsuda, S. Komiyama, T. Terashima, K. Shimura, and 
Y. Bando, Phys. Rev. Lett. 69, 3228 (1992). 

46 J. C. Garland and H. J. Lee, Phys. Rev. B 36, 3638 (1987). 

47 A. F. Hebard and A. T. Fiory, Phys. Rev. Lett. 50, 1603 
(1983). 

48 H. S. J. van der Zant, H. A. Zant, and H. A.Rijken, and 
J. E. Mooij, J. Low Temp. Phys. 79, 289 (1990). 

49 a) S. T. Herbert, J. Markus, S. E. Hebboul, and J. C. Gar- 



18 



land, J. Low Temp. Phys. 106, 607 (1997); b) S. T. Her- 
bert, Ph.D. thesis, Ohio State University, 1993. 

50 D. C. Harris, Ph.D. thesis, Ohio State University, 1989. 

51 D. W. Abraham, C. J. Lobb, M. Tinkham, and T. M. Klap- 
wijk, Phys. Rev. B 26, 5268 (1982). 

52 D. W. Abraham, Ph.D. thesis, Harvard University, 1983. 

53 J. Maps and R. B. Hallock, Phys. Rev. B 27, 5491 (1983). 

54 R. J. Dionne, S. B. Crary, and R. B. Hallock, Phys. Rev. 
B 33, 1997 (1987). 

55 S. W. Pierson, Phys. Rev. B 51, 6663 (1995). 

56 M. Friesen, Phys. Rev. B 53, R514 (1996). 

57 P. Minnhagen and P. Olsson, Phys. Rev. B 45, 10557 

(1992) . 

58 V. Cataudella and P. M. Minnhagen, Physica C 160, 442 
(1990). 

59 S. W. Pierson, Phys. Rev. Lett. 75, 4674 (1995). 

60 M. Friesen, Phys. Rev. B 51, 632 (1995). 

61 M. Friesen, Phys. Rev. B 51, 12786 (1995). 

62 S. Sujani, B. Chattopadhyay, and S. R. Shenoy, Phys. Rev. 
B 50, 16668 (1994). 

63 K. Epstein, A. M. Goldman, and A. M. Kadin, Phys. 
Rev. Lett. 47, 534 (1981). 

64 D. P. Norton and D. H. Lowndes, Phys. Rev. B 48, 6460 

(1993) . 

65 P. M. Minnhagen, Phys. Rev. B 23, 5745 (1981). 

66 D. T. Fuchs, E. Zeldov, M. Rappaport, T. Tamegai, S. Ooi, 
H. Shtrikmem, Nature 391, 373, (1998). 

67 Shaw et alx3 do extract a value of z by assuming that the 
resistance R(T) cx £~ 2 . (See paragraph around Eq. (8) in 
that reference.) This is tantamount to assuming z — 2. See 
Eq. ©. 

68 S. Kim and M. Y. Choi, Phys. Rev. B 48, 322 (1993). 

69 See, e.g., Fig. 29 of G. Blatter M. V.Feigel'man, 
V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Rev. 
Mod. Phys. 66, 1162 (1994). 

70 S. W. Pierson and O. T. Vails, unpublished. 

71 V. A. Gasparov, I. E. Batov, C. Kwon, and Q. Li, preprint. 

72 K. Holmlund and P. Minnhagen, Phys. Rev. B 54, 523 

(1996) . 

73 R. H. Koch, V. Foglietti, W. J. Gallagher, G. Koren, 
A. Gupta, and M. P. A. Fisher, Phys. Rev. Lett. 63, 1511 
(1989). 

74 H. K. Olsson, R. H. Koch, W. Eidelloth, and R. P. Rober- 
tazzi, Phys. Rev. Lett. 66, 2661 (1991). 

75 K. Moloni, M. Friesen, S. Li, V. Souw, P. Metcalf, L. Hou, 
and M. McElfresh, Phys. Rev. B 56, 14 784 (1997). 

76 K. H. Lee and D. Stroud, Phys. Rev. B 46, 5699 (1992). 

77 H. Weber and H. J. Jensen, Phys. Rev. Lett. 78, 2620 

(1997) . 

78 S. E. Hebboul, preprint. 



19 



